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I show how the Fermi and Bose pressures in quantum systems, identified in standard discus- 
sions through the use of thermodynamic analogies, can be derived directiy in terms of the fiow 
of momentum across a surface by using the quantum mechanical stress tensor. In this approach, 
analogous to classical kinetic theory, pressure is naturally defined locally, a point which is obvious in 
terms of the stress-tensor but is hidden in the usual thermodynamic approach. The two approaches 
are connected by an interesting application of boundary perturbation theory for quantum systems. 
The treatment leads to a simple interpretation of the pressure in Fermi and Bose systems in terms 
of the momentum flow encoded in the wave functions. I apply the methods to several problems, 
investigating the properties of quasi continuous systems, relations for Fermi and Bose pressures, 
shape-dependent effects and anisotropics, and the treatment of particles in external fields, and note 
several interesting problems for graduate courses in statistical mechanics that arise naturally in the 
(~| ' context of these examples. 
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I. INTRODUCTION 



C/5 ■ 

^ , The concept of pressure in quantum systems is usually introduced in equilibrium quantum statistical mechanics using 
thermodynamic analogies. Thus, the Helmholtz free energy F is identified with the logarithm of the the canonical 
partition function Z, and the thermodynamic potential 51, with the logarithm of the grand partition function Z 
■ through the relations 



F^-kT\nZ, n^-kTlnZ. (1.1) 
Here Z is defined as the usual sum over the energies in the system, 

Z= Tre-'^« = ^e-'^^° (1.2) 



^""^ ■ with (3 — 1/kT. In the case of particle systems, a separate partition function Zjv can be defined for each particle 
' number N . Z is then defined for an indefinite number of particles by 

Z = Y,e'"'^ZN, (1.3) 

^ ' N 

^ ' where /i is the chemical potential. The pressure is customarily determined through one of the standard thermodynianic 

• ^ ' relations 
C/2 ■ 

P = -dF/dV ^{kT/Z){dZ/dV) or PV ^ -Vl, (1.4) 

O ^, with F evaluated at fixed temperature T and particle number N , and il, at fixed T and /i. In the second case, the 
IL* ■ average particle number and chemical potential are related by the condition N — —dfl/dii. 
. !^ ] These relations can be checked in classical particle statistical mechanics by an appeal to the results of kinetic theory, 
■ and can be further motivated by an appeal to the concept of generalized forces when the energy of a system depends 
] explicitly on external parameters such as its volume. However, standard discussions do not show directly how the 
pressure relations arise in a kinetic theory-like context in quantum statistical mechanics. For a sampling of standard 
treatments, see Refs. [lll2l3ll4ii5i.6«7„8ii9t 

The objective of this paper is to give direct derivations of the pressure in Fermi and Bose systems using ideas 
analogous to those in kinetic theory, specifically the relation of pressure to momentum flow and the quantum stress 
tensor. These are well defined ideas in quantum systems, give a "quantum kinetic theory" approach to pressure, and 
lead to a direct interpretation of the pressure in Fermi and Bose systems in terms of the momentum flow encoded 
in the wave functions. I will show, in fact, that pressure is naturally defined locally, a point which is obvious in the 
stress-tensor approach but is hidden in the usual thermodynamic approach. The two approaches are connected by an 
interesting application of boundary perturbation theory for quantum systems. 

The basic ideas and relations in my approach are developed in Sec. ^ and their use for systems of noninteracting 
fermions and bosons, in Sec. lIIII I then consider several examples including the properties of quasi continuous systems, 
relations for Fermi and Bose pressures, and shape-dependent effects and anisotropics in Sec. II and examples for 
particles in external fields in Sec. IIVBI The results lead to several interesting problems for graduate courses in 
statistical mechanics. 
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II. PRESSURE IN STATISTICAL SYSTEMS 
A. Pressure and the stress tensor 

The pressure of a system at a point a; on a surface S can be defined as the rate of momentum flow dp/dt across a 
surface element dS = fidS at x, 

P{x,t) = ^{x,t)-h, (2.1) 

that is, in terms of the force per unit area or stress acting across the surface. P depends implicitly on the orientation 
of the surface through h, but I will not indicate this explicitly both for notational simplicity, and because the apparent 
n dependence is, in fact, absent for the familiar extensive, quasi continuous systems. I will apply this definition in the 
quantum context. For definiteness, consider a quantum system of identical particles with the Lagrangian density 

c^'- {rdti^ - dtr i')-7r- Ec^"'//) • (v„^) - r^i^, (2.2) 

2 2m ^ — ' 

ji=i 

where = i){xi, . . . ,Xf^,t) is the many-particle wave function for the system and V(a;i, . . . , Xj^) is the potential. 
The corresponding Schrodinger equation is 

^2 N 

ihdt^ = -—Y,'^li' + y^- (2.3) 

n—l 

The momentum density at a point x for the system in a state \a) with wave function ipa is 

Pa{x, t)^^Yl / (^a^"^" - '^"C V-a) S^{X - X^)d^Xi ■ ■ ■ d^ X N , (2.4) 

a result obtained by integrating the sum of single-particle momentum operators {h/2i) {ipa^n^^a ~ ^n^^ai^a) over 
the coordinates of the unobserved particles. The integrations are over the volume V in which the system is confined. 
Similarly, the local number density is 

N 
n=l 

After some rearrangements, the time derivative of Pa can be written as 

N 



N 

{x,t)^Y^ J ij:li>a5^{x^x^)d^xi---d^XN- (2.5) 



^{x,t) - / (9tCV„^c- V„V:5tVa)^'(a;-a;„)d=^xi---d3 

n ^ r 

+ YiY. (C^tV-a - dt^c S^'ix - X^)d'x, ■ ■ ■ d^XN. (2.6) 

n—l 

It is convenient at this point to switch to a component labeling of Pa and consider dtPa,i- Using the Schrodinger 
equation H2.3|l to eliminate the time derivatives in H2.6|l . splitting the double sums that appear into terms with identical 
and different particle labels, and organizing the results as much as possible into a set of divergences, I find after a 
straightforward calculation that 



dt 



XN 



{X, t) = ~flj^*c {'^n,^VHaS'{x - a;„) d'xi ■■■d 
N 3 . 

E E / V„,fc (V„,jVa^n,fc'0a + V„,fcV'* V„,i-0Q) {x - X„)d^Xi ---d 

1 J, 1 



2m 



n=l k=l 
-2 N 



^ 53 / ^"'^ (^".fcC V„,fcVa) S^{x - Xn)d^Xi ■ ■ ■ d^XN (2.7) 
n—l 

- N 

J (C V?>a + Va) S^{X ~ Xn)d^Xi ■ ■ ■ d^XN 



h 

4m 



n=l 

surface terms 
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The surface terms result from the integration of divergences V; • (•) in variables xi other than the selected variable 
Xn = X using Gauss' theorem. These terms vanish for the usual boundary condition for the energy eigenstates needed 
below, that ipa = for any of the coordinates on the boundary of the confining volume V , and will be dropped. 

The first term on the right hand side of H2.7|l is just the force density Fa^i{x,t) at x. The remaining terms are in 
the form of a divergence, and the result can be written as 



dt 

or, in dyadic notation, 



'x,t) ^ F^,,ix,t) + VkT^^^ix,t), (2.8) 



^{x,t)^F^{x,t) + V-T''ix,t) (2.9) 



where T^j is the quantum stress tensor evaluated in the state \a), 

dC dC 

d[dk^*) d{dk'ip) 

= X! / [^n.iV'a Vn^fcV^Q + Vn^fcV'Q Vn^iV'a - 4,iV„V* • V„V'q (2-10) 

n=l 

-\^k,^ (C VnV'a + Va) {x " Xn)d^Xi ■ ■ ■ d^XN- 

Upon integrating H2.8|l over a volume V' C V , one finds that the total momentum p in V' changes both because of 
the bulk action of the forces, and from the flow of momentum across the boundary surface S' — dV' , 

dpa 



or, again in dyadic notation. 



^ix, t) = J^^ F^,,{x, t) d^x + j^^ dSk T^^x, t) (2.11) 



^{x,t)^ I F^{x,t)d^x+ I dS-T°'{x,t). (2.12) 



Here dS ■ T is just the rate of momentum flow across the surface element dS = ndS into the volume V' , with fi the 
outward normal to the surface. From (|2.1() . the pressure at x is given by the momentum-flow per unit area out of V' . 
That is, 

Pa{x,t) = -n-T°'{x.t)-n, a; G S". (2.13) 

For equilibrium quantum statistical mechanics, the relevant states |a) are stationary states, that is, energy eigen- 
states, with Vq = 4'a{xi, ■ ■ ■ , a;Ar)e~*^°*/^. In this case, the explicit time dependence drops out in the expressions 
(|2.7|) - (|2.14|) . and tpa can be taken in these and following expressions as the spatial wave function ipaixi, . . . ,xn). 
The pressure, stress tensor, and force density are then independent of t, Pa{x,t) — > Pa{x), T"-{x,t) T"-{x), and 
Fa{x, t) Fa{x). Furthermore, dpa/dt = so —V T = F, and the divergence of the local stress is balanced by the 
force density. I will specialize to this case for the remainder of the paper and use the definition 

Pa{x) = -h-T"{x)-h, xeS'. (2.14) 

for the pressure at x in the state |a), where T°'{x) given by H2.11|l with ipa the spatial wave function. 

I will first consider the case in which T^' = ^ is the volume in which the system is confined, and will consider a 
more general case in Sec. IIIII For x on the boundary surface S = dV, ipa and the derivatives of ip parallel to the 
surface vanish. Thus, using (|2.11|) . 

Pa{x) = ^Y1 y • '^'V':)(" • Vi4>c.)S^{x - xi)d^xi ■ ■ ■ d^XN, (2.15) 
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and Pa{x) depends only on the normal derivatives of a.t x. Finally, weighting Pa by the statistical factor e '^^^ 
and averaging over all energy eigenstates |a), 

t2 ^ /• 

^^^^ = Z 2^ ^ ^ / ^" ■ '^'^"^(^ • '^iMe-"''- S'{x - xi)d'xi ■ ■ ■ d'xN, (2.16) 

™ a 1=1'' 

where Z is the canonical partition function in 1)1.2(1 . The sums are over all completely symmetric states for Bose 
systems, and over all completely antisymmetric states for Fermi systems. 

I would emphasize that the wave functions ip(x)i, . . . ,xn should include any factors such as spin eigenfunctions nec- 
essary to describe internal structure that does not affect the original Schrodinger equation. The sums over eigenstates 
include sums over the extra quantum numbers necessary to label the states completely. In the case of observables 
such as the total number density or pressure that do not depend on the internal structure, the ip^s can be reduced to 
spatial wave functions, and the right hand sides of (|2.16() multiplied by the appropriate degeneracy factor g. I will 
follow this convention throughout the paper. 

Eauation H2.16|l . or the more general form in ((2.14(1 . gives my basic "quantum kinetic theory" result. The differences 
between fermions and boson enter only through the symmetry properties of the wave functions and the resulting 
differences in the sums over states. Before going on to investigate these in simple cases, I would reemphasize that the 
kinetic definition of the pressure is intrinsically local and is expressed through the action of the momentum operators 
—ih'V as would be expected on the basis of classical kinetic theory. It is not immediately clear how this definition of 
the pressure is connected with the usual "thermodynamic" definition in 1(1. 4() . I will first show that the two definitions 
are equivalent when one considers local variations of the volume in the relation P = kT din Z/dV. 



B. Pressure from the partition function 

The thermodynamic definition 1(1.4(1 of pressure in terms of the canonical partition function gives the relation 

„ kTdZ 1 dEa aj, 

a 

A comparison of this expression with ((2.16(1 suggests that dEa/dV should be expressible for local variations in V 
in terms of the normal derivative of ipa on the boundary surface. This is, in fact, easy to show using boundary 
perturbation theory. I consider a small change in the volume of the system implemented by moving the boundary 
surface outward over a small surface patch AS* through a normal displacement Sx = nSx(x) that varies smoothly 
over dS and vanishes elsewhere. The energy E'^ of the system in the distorted volume will differ from the energy E^ 
of the original system, with E'^. — Ea + 6Ea ■ The perturbed spatial wave function ip^ and the original wave function 
Ipa satisfy the time- independent versions of the Schrodinger equations ((2.3(1 . 

1=1 1=1 

Multiplying the adjoint of the first equation on the right by and the second equation on the left by V'q* and 
subtracting gives the expression 

6Ea i^':^Pa = E • (^'V'L* - V/^a) , (2.19) 

2m -^-^ V / 

1=1 

where I have assumed that the potential V is unchanged. An integration of all the coordinates over the original 
volume gives 

SEa / V'lVad'xi • • • d^XN = -— E / ^*^(^) • / (^'^" ~ '^'^") '^'(^ ~ ■ • ■ (2.20) 

where the surface integration is defined in terms the variable x. 

The wave function ipa vanishes on the original surface S, so the first term in parentheses in ((2.20(1 vanishes. 
Similarly, tp* = on the distorted surface S' so the second term vanishes except on the patch AS where the two 
surfaces differ, and the surface integration reduces to the patch A5. For small normal displacements Sx = h6x{x), 
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can be approximated to first order on AS using the first nonzero term in its Taylor series expansion relative to S' , 
« —5x ■ Vi/'q ~ — (n • Vt/'a) Sx, where I note that the variation of ip'^ for small displacements parallel to AS does 

not contribute to first order. Finally, taking ^p'^ equal to ipa in leading order and using the normalization condition 

for the wave function, I find a first-order expression for SEa, 

^ r r 

6E^^-—y dS5x {n-Vii;D{n-Vii;D5\x-xi)(fixi---(fxN. (2.21) 
2™ ^ Jas Jv 

The integrand in the volume integral can again be taken as constant to leading order for x on AS* . The remaining 
surface integral simply gives the volume change 5V = J^g dS 5x, so 

BF h'^ ^ r 

" "2^ ^ Iv^ ■ ^'"^"^^^ ■ '^'^"^^'^"^ " ^'^'^'^^ ■ ■ ■ ^^-^^^ 

Substitution of this expression in H2.17|l reproduces H2.16|l . and the kinetic and thermodynamic definitions of the 
local pressure P{x) agree. For homogeneous isotropic systems, the factor {n - 'Vip^){fi ■ V-j/'q) in the integral in (|2.21|) 
has the same value at all points on S, the surface integral can be extended to the entire surface, and the calculation 
reproduces the usual a;-independent expression P = kT din Z/dV. 



C. Generalizations 



1. Electromagnetic interactions 



The results above can be generalized in various ways. For example, in the presence of electromagnetic interactions, 
the Lagrangian in H2.2|l becomes 



n=l 



(2.23) 



where A and $ are the vector and scalar potentials, e is the particle charge, and the particles are again treated as 
identical. The pressure of the system is related to the rate of change of the kinetic momentum "mw"= p—^A summed 
over the particles. A rather lengthy calculation gives the result 



n=l 

where F is the Lorentz force density 

TV 



^ d e ^ 



(2.24) 



Fo,{x,t) = J2 I eij*^Ei^aS^{x-x^)d^xi---d^XN (2.25) 

n=l 

N 

+ ^ J ^(j,hVn - -A^ Tpl X Blpa - X {iTlV - - A^ Va 5^ [x - Xn)d^Xi ■ ■ ■ d^ X N , 



and is the gauge-invariant stress tensor, 



^ k,i — 



n=l 

-5k,i (ihV 



(2.26) 



S^{x — Xn)d'^Xi ■ ■ ■ d^XN- 



The force density has the expected form, F ^ eE + {e/c)v x B. 
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The pressure is given by (|2.13() . or in the equihbrium case of static fields A and $ and stationary states |a), by 
(|2.14l) . If the latter is evaluated on the confining surface where the spatial wave function ip^ vanishes, the A-dependent 
terms in (|2.26(l all drop out, and the time- independent pressure Paix) is again given by the the expression in (|2.15(l 
and depends only on the normal derivatives of ip and ip*. The more general expression in (|2.14|l can be used for a 
surface element in the interior of the confining volume. 

A calculation similar to that in Sec. IIIBI also reproduces H2.15|l for static fields and stationary systems. The 
A-dependent terms again drop out on the confining surface, and the expression for 6Ea obtained in boundary per- 
turbation theory reduces to H2.21|l . 

2. Wave-type equations 

Quantized systems of bosons such as photons, mesons, or phonons in a solid satisfy wave-type equations, for 
example, the standard wave equation 

^d^^-Y^VlcP+'^cP, (2.27) 

n— 1 

with or without the extra mass term or potential m^{xi^ . . . ^x^). The normahzation of (j) is given in terms of the 
covariant current density/*^ and reduces for positive-energy eigenstates to 

- {(p*dt<l> ~ dtcj)* <l>) d^xi ■ ■ ■ (fxN = 4>*4>d^xi ■■■d^XN = 1. (2.28) 

An appropriate Lagrangian density for such systems is 

/ 1 

C = hci -dtr dtCP - ^n'P* ■ V„<^ - 

\ 71=1 

The corresponding momentum density p(a;, t) is 



(2.29) 



pix, t) = — / C^"-^* + 9t<P* V„0) S^{x - Xn)d^xi ■ ■ ■ d'^XN. (2.30) 

^ n=l •' 

Calculations of dp/dt similar to those above give the same formal result for the pressure as in (|2.14f) . but with the 
stress tensor now given by 

TV 

Tk,^=-hcY / iVn,,^'* Vn.k(l> + 'Vn.k(l)* Vn^^(l>- Sk^.C) S^{x-Xn)d^Xi---d^XN. (2.31) 

n=l 

This can be written in an energy eigenstate \a) with wave function (f> = (paixi, . . . , XM)e.~^^'^*'/^ as 

5^{x-Xr,)d^xi---d^XN. (2.32) 



n = l • 



2 

The corresponding pressure for a; on a boundary surface S — dV where (j) = is given by 

N 

Pa{x)^hcY / {n-Vn(l)*a){n-Vn(l)a)S^{x~Xn)d^Xi---d^XN, (2.33) 

n=l 

and depends only on the normal derivatives of on the boundary. For a general surface S' C the full form of Tj^^ 
must be used. The pressure in the canonical ensemble is just 

P(^) = l^Pa(a;)e-'^^°. (2.34) 
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It is easy to show that the same resuh follows from the usual thermodynamic relation in H2.17|l . Thus, from H2.27|l . 

TT^'t'a = -V H -Q— '/'a- (2.35) 

h n 

An equation of the same form with a perturbed energy E'^ holds for the perturbed wave function which results 
from a local displacement of the boundary. Combining the two equations, I find as the analog of H2.19|l that 

2Ea SEa I , , 



Tic 



ncY,'^i-(yi'p<:<i^c.-'i^:vi(^S)- (2.36) 



1=1 



Manipulations equivalent to those following (|2.19|) and the use of the normalization condition (|2.28|) then give the 
result 

BF ^ /■ 

" {X) = -hcJ2j (f^- ViC) (" • Vi0a) S^{X - Xi)d^Xi ■ ■ ■ d^XN (2.37) 



dV 



for local variations of the boundary surface. Finally, the use of the thermodynamic relation 1)2.17(1 reproduces the 
expression for the local pressure in (|2.34|) . I would emphasize, however, that the equation Pa{x) = —n -Ta ■ n for the 
local pressure in a state |a) holds more generally than (|2.33f) and does not require that x be on the boundary surface. 

III. PRESSURE IN NONINTERACTING BOSE AND FERMI SYSTEMS 

A. Bose systems 

As a first example, I will consider the important case of noninteracting bosons in an external field . Interparticle 
interactions can be treated approximately using standard methods. See, for example, Refs. IHI^I^I H^I^ r^l^ 

The Hamiltonian for N noninteracting subsystems or particles is a sum of N identical single-particle Hamiltonians 
Hi, H = ^fLiHi{xi). The wave functions ipk for the single-particle states \k) satisfy the Schrodinger equations 
Hiipk{xi) = Ek^k{xi)- I will suppose that the energy eigenvalues have been ordered so that Ei < E2 < E3 < . . .. 
The total energies are simply sums of single-particle energies Ek , and can be labeled by the number of particles in 
each single-particle eigenstate 

Eni,n2,... ^ niEi + n2E2 + ■ ■ ■ whcrc ni+n2 + --- = N. (3.1) 

The full wave function for N bosons with ni in state |fci), n2 in state |fc2), ... is then a fully symmetric sum of product 
wave functions, 

111 factors n2 factors 

V'«i,n2,...(a;i, • ■ • ,a3Ar) = ^= X! [ V'l ' ' ' '/'l "02 ' ' ' "02 •••] {P{xi, . . . , Xn)) , (3.2) 

where the sum is over all permutations P of N objects, and wave function factors with nk — are to be replaced by 1. 
The coordinates of the successive wave functions with Uk ^ are given in each term in the sum by the corresponding 
coordinates in the permutation P of indicated. The set of n^'s gives a unique labeling of the state. 

Their values are restricted by the condition — N . 

The number density of particles at a point x is given for a definite state by H2.5|) . The total density reduces after 
the integrations to 



nnun2,...{x) = ^nkil)l{x)ipu{x), (3.3) 

k 



and a final integration over x gives the total number of particles or independent subsystems N since "Y^^^ Uk — N . 
Similarly, from H2.15|l . the pressure on the boundary surface associated with the given state is 

Pni ,n2 , . . . ("^ 

) = ^nfe(n • Vrk{x)){n ■ V^Jkix)). (3.4) 
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More generally, for a surface S" — dV , V C V, 

Pn,,n-,,...{x) = -^nkn-T'' -n, x e S' , (3.5) 

k 

where the tensor Tj^j is given in H2.11|) . 

It is difficult to work with the canonical distribution for bosons because of the restriction J2k n-k ~ N. I will 
therefore change to the grand distribution as is usually done. Multiplying H3.3|) by the Boltzmann factor e~'^^"i'"2'' ■ 
for the specified energy and by a factor e''^^ which will be used to enforce the correct average number of particles, 
and summing over the and A^, I find that 



z 

k ni,n2,... N 

^E (E'^'=^fe(^)^fe(^) 

k \ rik 



-f3{Ek~tj.))2 1 
k ^ ' k'^k 




where convegence of the sum requires that Ek — > 0. Z is the grand partition function. 



Z 

k 



so the result reduces to 



n{x)=Y,rk{x)Mx){e'''^'''-'''^ (3.8) 



This is just the result that would be expected. Each product of wave functions V'fcV'fc appears with a weight which 
is just the average occupat ion n umber of the state \k) as calculated for the usual Bose distribution for noninteracting 
particles as in Refs. llOalalil lsla lTllsllgi 

In the presence of spin or other internal degeneracies, the ?/;'s can be reduced to the spatial factors in the full wave 
functions, and after multiplying the right hand side of the equation by the appropriate degeneracy factor the sum in 
(|3.8|) can be taken to run only over nondegenerate energies with the internal factors in the wave functions suppressed. 

The integral of n{x) gives the average number of particles N in the entire distribution, 

7V = ^(e'^(^--^)-l)"\ (3.9) 
fc 

This gives an implicit relation for ^ in agreement with the thermodynamic expression A'^ = kT{dhiZ /dfi). For T ^ 0, 
^ ^ Eo — ^ + ■ • • , and the particles collect in the ground state |0) with n{x) ~ N'i(jq{x)iPo{x). 

A similar calculation gives the result for the pressure at x, again of a form that could be anticipated from the 
single-particle forms of H2.11|l . (|2.14ll . and (|2.15|l . 

P{x) = -^n-r'^ •n(e'^(^^-^) - l)"^ (3.10) 

fc 

= 1^ E • ^^*k{x)) • ^M^)) {e^^''--'^^ - l) . (3.11) 
fc' 

The first form holds for a; on a surface S' inside the confining volume, the second form for x on the confining surface 
S = dV. For 0, (Tmil gives P{x) {NTi^ /2m){h ■ V %{x)) {h ■ V ^Ijo{x)) . 

It is interesting to note that the pressure does not vanish exactly for a system confined in a finite volume even at 
T = 0, a result connected to the kinetic picture and the uncertainty relation. For example, for Bose particles in a one di- 
mensional box of length L, the wave function are -i/'n = \/2/Lshi {nnx/L), and H2.15|l gives P„ — {N/L){h'^TT^n'^ /rriL^), 
n — 1,2, . . .. All the particles collect in the ground state at T = (Bose condensation), but the pressure is nonzero 
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for finite L, vanishing as L^^ for L — s- cx) at fixed particle density Nj L just as would be expected from the classical 
kinetic picture and the uncertainty relation p oc T%/ L. 

The pressure in systems of noninteracting bosons satisfying the wave equation (|2.27(l is also given formally by the 
expression in H3.1U|I . with T now given by the single-particle form of H2.33|l . The analog of H3.11|l is therefore 

P{x) = hcY^ (h ■ V0^(a;)) (h ■ VMx)) (e''(^'=-'^' ^ ' (^'^^^ 
fc 



B. Fermi systems 



The composite state of N noninteracting fermions is specified completely by giving the number of particles in 
each completely labelled single-particle state \k) where = or 1 only. The energy of the state \ni, n2, ■ ■ ■) is just 
En-i.n2,... — X]fc"-'=^'c in l|3.1|l . The corresponding wave function is given by the completely antisymmetric sum 

III »12 

V'«i,n2,...(a;i, • ■ • = ^=X!(^)^[ ^1 ^2 ■■■]{Pixi,...,XN)), ^nk=N, (3.13) 

where (— )^ is the signature of the permutation P of N objects. The factors of wavefunctions -0^ with = are to 
be replaced by 1. The coordinates of successive wave functions with Uk — 1 are given in each term in the sum by 
the corresponding coordinates in the permutation P of Xi, . . . ,xj^. The n's are restricted by the condition indicated, 
that their sum be N. 

The number density and pressure of the particles in the specified state are given at a point x by H2.5|l and H2.15|l . 
respectively, and reduce after the integrations are performed to the expressions in (|3.3() and (|3.4() . or more generally 
(|3.5|l . just as in the bosonic case. The difference between the two cases is entirely in the allowed values of the n's. 
The fermionic sum can be performed simply in the grand statistical distribution, with, for example. 



fc ni,n2,---=0,l N 

fc k'^k 

= ^ rkix)Mx) (e''(^'=-^) + iy\ (3.14) 



fc 

where I have used the relation, also easily derived. 



Z = [](l + e"''(^'="'^)) . (3.15) 
fc 

The expression in H3.14|l is again what would be expected since the final factor is just the average occupation number 
of the state \k) in the grand ensemble. Similarly, using H3.4() . 

^(^) = 1^ E • ^^kix)) (n ■ V^{x)) (e'3(^-Ai) + -\ (3.16) 
fc 



or more generally. 



P(a;) = -^n-T'^ •n(e^(^^'-'^) + l) \ (3.17) 



where T*^ is the single-particle version of (|2.11|) . In the case a system has internal spin-type degeneracies, the ip^s 
in (|3.17|) can be reduced the spatial factors in the full wave functions, and the sum restricted to the nondegenerate 
spatial eigenstates after multiplying the right hand side of the equation by the degeneracy factor g, 
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IV. EXAMPLES 
A. Quasi continuous systems 

1. General constderations 

The systems to which statistical descriptions are apphed most frequently are large, extensive systems in which the 
potentials are uniform or periodic. A well-known theorem shows that the number of eigenvalues Ek smaller than a 
fixed value E grows proportionally to the volume of the system for ^ — > oo. [See Kac, Ref. fill for a famous discussion 
of this result and its history in the context of the spectrum of a drum.] The eigenvalues therefore pack together for V 
large, surface effects on the spectrum become negligible, and it is plausible that the sums over states in the preceding 
sections can be converted to integrals when there are many states with energies less than kT . The main question 
concerns the behavior of the wave functions in the limit of large V . We expect, in fact, that the products i})*{x)iIj{x) 
and W%Ij*{x) ■ 'Vip{x) will each reduce for large V to the sum of a term describing their smooth average behavior, and 
extra rapidly oscillating terms that average approximately to zero. The result should again be insensitive to surface 
effects for V sufficiently large. 

These ideas can be illustrated for a uniform system in D dimensions by quantizing in a box with sides Li, i — 
1,2, . . . , D. The wave functions and energies are 

I 2 ITHj ' X ^ Jlj^ 

^{n} (^) - n V I" ' ^{"> -T.<^^ (4-1) 

where {n\ is a multi index, — (ni, . . . , njj) with — 1, 2, . . . only. This gives 

^|„^(a.)V.M(^) = (l-2cos^^) , (4.2) 

i—l ^ * / 

where Vd = III volume of the D-dimensional parallelepiped in which the system is confined and the state 

label k in earlier equations is now given explicitly by the multi index ni, 712, . . . , nu. 

The typical index for states excited at temperature T is n^x ~ {SmkTV'^/ ^ /K^Y/"^ . If this is large, many states 
will be excited as required for the conversion of sums to integrals, and the oscillating terms in (|4.2() will average to 
zero over small regions of the box. Then, for observations over such regions, V'*^ ~ l/l^ni ^ result independent of 
the shape of Vq. This is the same result as that obtained using the standard approximation of running waves with 
periodic boundary conditions, ~ {1/^/V) exp (^^ TrriiXi/ Li). An independent argument shows that the sums of the 
oscillating terms vanish rapidly at fixed x as the numbers of significant terms in the summations grow, that is, for 
many states excited. Since the level spacings tend to zero for Vu 00, either argument shows that only the leading 
term in 14. 2|) is important for spatially large systems. 

Dropping the oscillating terms in (|4.2(l . the expression for the local number density for uniform Fermi and Bose 
systems becomes 

E (e^(^"- --^^il)"', (4.3) 



where the upper and lower signs refer to Fermi and Bose systems, respectively. The sums can be converted approxi- 
mately to integrals by repeated use of the Euler-Maclaurin summation formula 



J2nn) = J^ fi^)dn-^f{Q)-^nQ) + j:^f (0) + - 



(4.4) 



The odd-order derivatives that appear in the Eulcr-Maclauring formula all vanish at rij = 0, 00 for the function 

in (|4.3(l . The first two terms in l|4.4|l are therefore all that survive up to exponentially small corrections that can be 
investigated using Poisson summation. Retaining only the leading corrections. 



n[x 



\- dn.dn, ...dnn (^^(S.C'Vs^L?)"?-.) ± l) "7l - i ^ <5(n,) + ■■■). (4.5) 
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At this point, a change to the momentum variables pi — {h/2Li)ni gives the famihar expression for the leading 
term, plus corrections that vanish as Vjj for Vd — > oo, 

nix) = / ^ (e/5(py2™-^) ^^y'_Y^±_l (^^Pipl_J2m-^.) ± l) + 0(1/ L^). (4.6) 

The momentum integrations extend over the infinite interval (— oo, oo), a Brillouin zone, or otherwise as appropriate. 
The corrections are of order h/ptLi for pi the typical value of the component of the momentum in the leading 
term, that is, of order Ti / \/ mkT Li for nondegenerate systems. 

Note that the final result for the leading term is isotropic in momentum space even though the original spectrum in 
(14. 1|) is different for motions in the different directions. This appears to be general for quasi continuous systems; see, 
for example, Sec. lIVBTI With enough energy levels occupied in the thermal distribution, the details of the spectrum 
become unimportant. 



2. Fermi and Base pressures 

A calculation of the pressure using the method above and either of Eqs. (|3.11(l or H3.12|l for Bose systems, or (|3.16|l 
or (|3.17|) for Fermi systems leads to analogous results for the pressure on a surface with normal h, 

J m \ / 

^ I / 7?^^^) (e/^t^-/^-^) ± l) ... = 1(B) + ... , (4.8) 

with corrections that again vanish as Vjj^^^ for Vd — > oo. The Bose and Fermi statistical factors in the integrals 
are isotropic in momentum space. Thus, the leading term in the expression for the pressure is independent of the 
direction of n, and the results can be expressed in terms of the average energies as indicated. 

I would emphasize that this result for the local pressure follows directly from the definition of the pressure in terms 
of the stress on a surface. The factors of n ■ p in H4.7I) arise from the momentum operators ~ih^ in —n ■ T ■ n, and 
correspond directly to the momenta that appear in the elementary classical derivation of the pressure in a gas. That 
is, the pressure is associated with the "beating of the particles against the wall." 

The thermodynamic definition gives the same pressure for the quasi-homogeneous system under consideration, and 
the two definitions are connected by boundary perturbation theory as shown in Sec. Ill Bl The derivation given there 
can be generalized to an arbitrary surface inside the volume V, but is only useful provided that, as here, enough is 
known about the wave functions to allow explicit evaluation of their derivatives. 

It is easy to derive the nonclassical properties of the Fermi and Bose pressures. An integration by parts in spherical 
coordinates brings 1)4. 8|l to the form 

In (l ± e-'^tPo-'^)) = TkT / In (1 T n{p)) (4-9) 

where n{p) is the Fermi or Bose statistical factor in (|4.8|l . Using the inequalities 

-ln(l-a;)>a:: and ln(l+x)<x (4.10) 

and the fact that the integral of n{p) gives the number density n{x), one finds that 

^'Formi(3;) > n{x)kT and FBoso(a^) < n{x)kT. (4.11) 

The difference clearly arises in the momentum-flow or stress picture from the necessity that the occupied single- 
particle states all be different for Fermi-Dirac statistics. This forces the appearance of higher momentum states than 
are needed in the Bose-Einstein case, and a higher pressure for fixed N and T. 

As an example of Bose pressure, I will calculate the pressure of an equilibrium system of noninteracting neutral 
mesons with mass m. The system will be taken as extensive or quasi continuous in dimension D. The pressure on 
the boundary surface is given in a state \a) by (|2.33() . It can also be calculated on an interior surface using the stress 
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tensor in H2.32|l and the definition in (|2.14|l . giving the same average resuh for large volumes Vd- Using (|2.33|l and 
the single-particle wave functions in a box normalized according to 




P.....P.v^) = ^^11^1- sin P^ = ^n., n. = l,2,..., (4.12) 

and averaging with the Bose statistical factor, I obtain 



where E{p) = y^p'^c^ + m^d^. In the limit to ^ 0, this reduces to the expression for the pressure for black body 
radiation or for phonons in a solid up to the necessary inclusion of the statistical factors for spins or polarizations and 
the use of the correct ranges of integration in the case of phonons. Thus, for black body radiation in three dimensions, 
including the spin degeneracy factor 2, and using the isotropy in momentum space and the fact that /i = because 
photon number is not conserved, 

the usual result. However, I would again emphasize that the calculation is direct, and is clearly connected to the flow 
of momentum across the surface through the discussion in Sec. Ill C II No thermodynamic relations were used. 

An example for Fermi systems that makes a good homework problem and shows the relation of pressure to mo- 
mentum flow is the calculation of the Fermi pressure at T = staring from the stress tensor. The usual argument for 
completely degenerate Fermi systems shows that all energy levels up to a Fermi energy Ep determined by N must be 
occupied. Then from l|TT7jl . P{x) = - Y,k d{EF - Ek)n-T^ -n where is the step function, 9{x) — 1 (0) for x > (<)0. 
The remaining calculation is simple for continuous systems, and leads directly to the thermodynamic result without 
the use of any thermodynamic relations. 



3. Anisotropic pressures 

The corrections from the conversion of sums to integrals in H4.6|l . and the corrections from finite-size effects in the 
spectrumjii are shape-dependent. This shape dependence leads for finite systems to anisotropic stresses or pressures. 
Consider, for example, the limit in which one of the dimensions of the rectangular box considered above, say Li, 
becomes small while the other dimensions remain large. If the lowest energy for motions in the 1 direction is large on 
the scale of kT, /SraL\ ^ kT, the sum over ni converges rapidly, and conversion of that sum to an integral with 
only small residual corrections is not possible. Keeping just the leading term in ni and treating the large dimensions 
in the continuum limit, the number density becomes 

Li Li J 11 \ / 

where /i' = /i — Ei with Ei =^ /%mL\ the ground state energy for motion in direction 1. Higher terms in the sum 
on ni are nominally suppressed by powers of exp[— (n^ — \)Ei/kT] ^ 1, ni > 1, but the situation becomes more 
complicated and some excitation must occur for Fermi systems in which the Fermi energy for N particles m D — \ 
dimensions exceeds Ei. I will not consider this refinement. 

The leading factor in 14.1511 is the absolute square of the normalized wave function 'ipi{xi) and integrates to unity. 
The momentum integral is independent of the remaining coordinates X2, ■ ■ ■ ,Xj), so an integration of n{x) over the 
full volume Vd gives the total particle number N as 

TV = Vd-1 J (^e-^ipl-i/^---^') ±iy\ (4.16) 



The number density can therefore be written as 
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The pressure on the wall of the box at xi = can be calculated using the expressions in lITT^ and lITTTIl . with the 
result 

N N 

Pi = -^-^^2Ei^. (4.18) 

All particles must be in the ni — 1 state in xi for kT <C Ei, with the effects of Fermi or Bose statistics absorbed 
in the integral factor in H4.15|l and the corresponding factor in the expression for Pi. The same result for Pi holds 
on the surface Xi = Li. Note that Pi is independent of the coordinates X2t ■ ■ ,xo that specify the location on the 
surfaces at xi = 0, Li on which the pressure is observed. 
In contrast, the pressure on any of the remaining walls is 

P. (-) = f sin^ ^ / L-^iA-.n-^'^ ± l) " - f ^ • ^f^lT^' (4-19) 

Li Li J h^^^ m \ J Li Li D - 1 Vd-i 

J > 1, where (E) d^i the average single-particle excitation energy for a continuous system in _D — 1 dimensions. The 
pressure is independent of the coordinates X2, ■ ■ ■ ,xd, but is modulated with respect to Xi by the factor |'i/;i(xi)p ~ 
(2/Li) sin^(7rxi/Li) which specifics how the N particles are distributed with respect to xi. The last factor in (|4.19|) 
is just the pressure in Z) — 1 dimensions. When the distribution in xi is not observed, xi can be integrated out. The 
integration reduces Pn to Pd~i, and the effective dimensionality of the system is reduced by one for kT ^ Ei, the 
limit in which the no thermal excitations in the 1 direction are possible. 

The effects are small in practice, with Ei/k = 0.2 K for a helium atom confined in a gap with Li ~ 1 nm. 
Anisotropics in the pressure would only be observable at lower temperatures. The situation is more complicated for 
electrons because of the effects of Fermi-Dirac statistics. 



B. Pressure and number density in an external field 



1. WKB approximation 

It is simple to treat the problem of otherwise noninteracting particles in a one-dimensional external potential V{x) 
using the standard WKB approximation discussed in most texts on quantum mechanics (see, for example, Ref. Il2i') . 
The single-particle energies iJ^ are determined in this approximation by the semiclassical quantization condition, that 
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k + -jh = 2j p{E,x)dx (4.20) 

for E = Ek, fc = 0, 1, 2, . . .. Here p{E, x) = ^2m\E — Y{x)\ is the local momentum defined classsically at energy E', 
h is Planck's constant, and xi, xi are turning points in the classical motion where p{E,x) = 0. The approximation 
can be shown to be good when there are many local wavelengths h/p between the turning points, but tends to be 
good even for low-lying states in the spectrum. Since the typical excitation energy in statistical systems is i? ~ kT, 
the WKB approximation will be valid provided J p(kT,x) dx ^ 1. Treating k and E as continuous, the density 
of states dk/dE implied by H4.20|l is 



dE hJ^VE- V{x 

The WKB wave functions can be written between the turning points asi^ 

o \ 1/4 



dk I 2m , , , 

dx. (4.21) 



M^) - ^-^—^j cos ^Cfc(x) - - j dx, (4.22) 
1 



^k{x) ^ - J y/2m[Ek - Vix)] dx, (4.23) 

and decrease exponentially outside that region. Ignoring the small contributions from the exponential regions and 
replacing the square of the cosine by its average value of 1/2 for many oscillations in the region in which E — V 
changes significantly, I find that the normalization constant is related to the density of states by 
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Thus, following the discussion of quasi continuous systems above, the local number density in the system is 




(4.25) 



where in the last two lines I have first replaced the square of the cosine in V'^V'fc by its average value 1/2 and converted 
the sum over k to an integral over E using 14.24|l . and then converted from E to p as the integration variable with 
p defined by the relation p^ /2m = E — V{x). The replacement cos^ ^1/2 may be taken as a local averaging when 
there are many oscillations in the region observed. Alternatively, I note that the zeros of successive eigenfunctions 
interweave, so the zeros in the individual terms in n{x) are washed out in the sum when many states are excited. 

The result in H4.25|l is just that obtained through thermodynamic arguments^ by dividing the system into small 
volumes over which V{x) can be taken as constant, and then considering the equilibrium of the subsystems. The 
result, as here, is to replace the chemical potential /i in the corresponding expression for free particles by /i — V{x). 
However, it is clear from the WKB-based derivation above that there are two key points in the quantum treatment. 
First, the average number k of excited states must be large enough and vary smoothly enough for energies on the 
scale of kT that the sum over states can be replaced by an integral. Second, the square of the wave function must 
oscillate sufficiently rapidly over regions in which E — V{x) changes significantly that the replacement cos^ ^ 1/2 is 
valid in the sum in H4.25|) . The second requirement is closely linked to the conditions needed for the validity of the 
WKB approximation, and for the replacement of the sum by an integral. It is worth noting in this connection that 
the wave function is defined over the entire volume in which the system is confined, and not just subvolumes as in 
Ref.B 

The expression in 1)4.25(1 can be extended immediately to three dimensions for systems with with additive potentials, 
V — Vi{xi). It can be extended to general potentials V{x) in the form 

n{x) - / § (e/?[pV2™-M+^/(x)] ± 1^ (4 26) 

using thermodynamic arguments as in Ref. 6., or directly using functional integral methods such as those in in Refs.lT^ 
andllJ. It is again required that the oscillations in if) be rapid on the scale at which E — V{x) changes significantly. 
The general result in 14.26|1 could probably also be derived in three dimensions using a WKB-like phase-integral 
approximation such as that investigated by Gutzwiller, Ref. Il5l but I have not attempted this. 

Finally, the total number of particles in the system is given by the spatial integral of n{x) over the confining volume. 



TV = 1 (^^P[v^/2.^-,+vi.)] ^ (4 27) 

giving a formula that can be used to determine /i. 

The pressure or stress in the external potential can be treated similarly. I will work in the interior of the total 
volume and use the general expression for the pressure given in ((2.14(1 , and the single-particle form of the stress tensor 
in 1(2.11(1 . specialized to one dimension. The key step involves the recognition that the derivatives in 1(2.11(1 can be 
taken to act only on the cosine factor in 7/;^, ((4.2211 . This gives 

^ - (^^) v/M^^^VMsin - 5) dx. (4.28) 

The term omitted is of relative order \^/%'k{E — V), and can be neglected in the region in which the WKB 
approximation is validji^ namely that the change in the potential over a wavelength X = h/p is small on the scale of 
E — V. Furthermore, the correction term oscillates out of phase with the main tern, and interference effects can be 
neglected in averaging V'^V'fc A similar result holds for the second derivatives, with the neglected terms just those by 
which the WKB wave function fails to satisfy the exact Schrodinger equation. 
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The result of the calculation is 



P{x) » ^AA;,V2™[Sfc - V{x)] (e/3(^^-^(-) ± 

dE^2m[Ek - V{x)] (^eP{E,-v(x) ^ (4 39) 

h m \ J 



The expression in H4.29|l be generalized to more dimensions using thermodynamic arguments or functional integral 
methods, and should properly be stated in terms of the stress across a surface with normal n, 



-fi - T ■ n 



m 



(e'5[p'/2"-'^+^(") ± 1) \ (4.30) 



The final factor in (|4.3U|) is isotropic in p so the stress at a given point x is the same in all directions even for V{x) 
anisotropic. Thus, (fi • p)^) = {p /D), and P{x) = j^K{x) where K{x) is the average kinetic energy density at x 
with the average taken over the local statistical distribution. 

I turn next to two examples which illustrate the effects of external fields in interesting physical situations. Both 
lead to useful homework problems for graduate courses in statistical physics. 



2. Example: Particles in a linear potential 

For particles in a linear potential V{z) = Vq x (z/zq) with no potentials for the motion in the transverse directions, 
the motion in z can be described in the WKB approximation as above, while the motion in the transverse coordinates 
can be described in terms of running waves with momenta pj^ = {px,Py)- The single-particle energies are 



2m zo \32mVoJ \ 4 

where the second term is the WKB energy of the vertical motion. Many states of the vertical moton will be excited at 
the thermal energy kT for {kT/Vo)'^^'^(8mVoZQ/9Tr'^h'^)^^'^ ^ 1, a condition always satisfied under realistic conditions 
for gases in a gravitational field with Vq — mgzQ or electrons in a constant electric field Eq, Vq = cEoZq. The sums 
over k can be replaced by integrals over a momentum pz defined to reduce the second term in H4.31|l to standard form. 



VJl - — _ _ 

4 J ' 2m Zo \32mVQ 

The corrections for Fermi or Bose statistics are unimportant for gases in a gravitational field under normal con- 
ditions. Following the development in Sec. IIVBTI with the Fermi or Bose factors replaced by the simple Boltzmann 
factor e^'^'^^'="''\ one obtains the classical barometric equations 

F(z) = F(0)e-™»"/'='^, n{z) = n(0)e-™»"/'='^, P(0) = n{0)kT. (4.33) 

For a system with area A, the number density n{0) at z = is given in terms of the total number of particles N by 
the integral 

d^xnM = n(0) / d^xe-'^a^^kT ^ ^.qw ^ ^4 34) 

Jo 

so rt(0) = (N/A){mg/kT) and P(0) — Nmg/A. The result is as expected. Note, as remarked above, that the pressure 
is isotropic at any point, Tx^x = Ty^y = Tz^z, even though the potential is not, and isotropy was not used in the 
derivation. This is general for quasi continuous systems at sufficiently high excitation. 

A more interesting result with respect to the gravitational field is the existence of a Bose-Einstein condensate in an 
ideal system at sufficiently low temperatures, an example that makes a good homework problem in a graduate course. 
The particle number for the Bose system is given in H3.9|) . This becomes 

3A /" ,9 f , 9 / rtr»,2 _ ,,] ^\ 1 iA 



pIk^ 



J d'p^J^ dpzpl (e'3[pV2™-.l _iy\^J d^p^J^ dpzpl [e^^'^"' - l)~' (4.35) 
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for the energy spectrum in (|4.31|l . with sums converted to integrals and Pz defined through (|4.32() . The inequahty 
follows from the convergence requirement that Eq — jj, > and the approximation in H4.35|l that i?o « 0. 

The inequality is clearly violated for a fixed N at sufficiently low temperatures or large values of It is then 
necessary to single out the ground state as this is given zero weight in the transition from a sum over states to an 
integral over E or p, and include its occupation number Nq explicitly. N is then given by 



N = No + 



3A 



d p± dpz p', 



/2mkT 



(4.36) 



The integrals can be evaluated by changing to spherical coordinates, and then to the variable t — p^ /2m after 
performing the angular integration. The final integral gives a product of a generalized factorial or gamma function 
with a Riemann zeta function , 



dt 



1 



r(.)c(^). 



(4.37) 



The inequality in (|4.35|) can just be satisfied for a given iV at a temperature T = Tc determined by setting the right 
hand term equal to iV, and fails at lower temperatures. The calculation gives 



N 



corresponding for a gravitational potential to 



Ap^h^ 



(2mfcTe)^/2C (5/2) 



mg / 
C(5/2) \2^J 



-1 2/5 



At lower temperatures, 



No = N 



1 - {T/T, 



n5/2 



(4.38) 



(4.39) 



(4.40) 



The power of T/Tc in (|4.40() is different from that for an ideal system with no field present. The critical temperature 
Tc is also higher for fixed N than the critical temperature in the absence of the gravitational field, Tc/T^ — 
[C(3/2)/C(5/2)]^/^ « 1.56, a difference attributable to the greater density of the gas near the ground. Finally, the 
ground state wave function is compact in z with a characteristic extent Zmax ~ Ei/mg = {81h^ /512mgY/^ w 5.4/im 
for helium, and the condensate "falls to the fioor." 

It is also interesting to note that the presence of a gravitational potential leads to the appearance of Bose-Einstein 
condensation in a two dimensional system with V = Vo{z/zo) and free motion in a box of length L in the transverse 
direction, with 



No=[l- {T/TcY 



kTr 



1 f2plhN 
m \ TT^ L 



1/2 



(4.41) 



There is no condensation for free motion in two dimensions. 

The calculations of the local number density and pressure in an external field are simple for dilute systems for which 
the Fermi and Bose statistical factors reduce in first approximation to the usual Boltzmann factor. They cannot be 
done exactly when quantum corrections are important, with many particles within a volume of a thermal wavelength 
cubed, but are similar numerically to the calculation of N{x) in the following example. 



3. Example: Bosons in a harmonic trap 



An example of interest in connection with Bose-Einstein condensation is that of atoms confined in a harmonic trap 
with V{x) — J2^=i • The total number of particles in the system is given by (|4.27f) . specifically. 



N = 



dPx d^p 



|^e'^i:jP?/2ni+(nic^?/2)a;?]-/5M _ \ (4 42) 
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The integral can be simplified by changing to the dimensionless variables x[ — LOi^ m/2kTxi, p'^ = pi/V2rnkT, and 
then going to six- dimensional coordinates s — {x' ,p'), — J2iipi^ + ^i^)^ ^^'^ working in a spherical representation. 
The result is 

N--^ ^] — . / dn,. (4.43) 



Ti J uJiLU2'^3 Jo e^^^f^f^ — 1 

Here dfig is the element of solid angle in six dimensions, and / dVL^ = tt"^. 

The remaining integral is maximized for /i = 0, and can be evaluated exactly in this limit by changing the integration 
variable from siot = s^ and using (|4.37() . The resulting equation determines the critical temperature for the onset 
of Bose-Einstein condensation for fixed particle number N , 



kT, = n \^ui^u;2CJ3^ j . (4.44) 

For a spherical trap with an oscillation frequenc y — 150 Hz and = 4 x 10'*, fairly typical conditions for 
original experiments with Rb atoms, Refs. 'l6*17'l^ this equation gives Tc = 6.77 nK x N^^^ = 232 nK. Note that 
tkij/k = 7.20 nK <^ Tc, so a large number of oscillator states are excited at Tc, and the use of the integral approximation 
to the sum over states is legitimate. 

The number of particles in excited states for T < Tc, ^ = is A^cxcitod = N(T/Tc)^, and the number in the ground 
state is therefore 

iVo = iV [1 - (T/Tcf] , T<Tc, (4.45) 

with Nq = for T > Tc. These calculations illuminate the conditions under which a real Bose-Einstein condensate 
can be formed in a gas, and make good homework problems. 
The number density in a spherical trap follows from (|4.42|) . 

n(r) = iVo(T)|Vo(r)P + ^ J d^^P (e^fp'/^""''-'/^) - l) (4.46) 
for T < Tc, where r — \fx^ and 

^o(^) = ^T7TTe~'''/'^o^ ro = V^/ww. (4.47) 

The actual evaluation of the local number density in the trap from (|4.46|l requires some numerical calculation but 
gives a striking illustration of the emergence of the condensate. It is useful to scale r by tq, n{r) = d^N/d^r by N, 
and T by Tc, and change to t = p^ /2m as the variable in the final momentum integration. Then, using (|4.44|l . H4.45|l . 
and l|TT7|l . 



1 d^N _ 1 
N'dHrJ^ ~ ^ 



To' ^ ■ 

1 f hujV'^ fT^ 



V2^2^(3) \kTcJ \T, 



dtVt(e*+^(''"/'=^^)(^=/^)('^'/''«)-l) (4.48) 



A very sharp ground-state peak appears in the initially rather broad in the number density as as T is lowered below 
Tc in qualitative agreement with the original experiments in This makes an an interesting comparison, and 

gives a real feeling for how the theory relates to observed Bose-Einstein condensates. 

The pressure in the trapped system can be calculated similarly, and balances the force from the confining oscillator 
potentials. However, because of the long mean free path for particle intractions, it is not relevant for the expansion 
of the condensate when the confining interactions are suddenly removed. 



V. COMMENTS 



The main objective of this paper was to give direct derivations of the pressures in Fermi and Bose systems using 
the relation of pressure to momentum fiow and the quantum stress tensor. This "quantum kinetic theory" approach 
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is simple conceptually, and shows that the pressure is naturally defined locally, a point of interest for particles in 
external fields. It leads also to a direct understanding of the difference in Fermi and Bose pressures at fixed particle 
number and temperature in terms of the different momentum states excited, a point often argued qualitatively. A 
bonus of the analysis was the appearance of the simple examples of the use the use of boundary perturbation theory in 
quantum mechanics necessary to establish the connection of the of the usual thermodynamic arguments for particles 
and fields to the results obtained directly in stress-tensor approach. 

I also discussed the properties of extensive, quasi-continuous systems, showed the role of excitations high on the 
scale of kT in obtaining isotropic pressures in intrinsically anisotropic systems, and illustrated the appearance of 
anisotropies and the effective reduction of the dimension of a system at low enough temperatures. Finally, I gave an 
explicit WKB derivation of the usual expressions for the number density and pressure of particles in an external field, 
and presented several examples which illustrate the use of the stress-tensor method in real physical problems. I have 
found these examples to make good homework problems in a graduate course on statistical mechanics. 
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